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Abstract 

We demonstrate that Plebahski's first heavenly equation decouples in infinitely 
many ways into a triple of commuting (1 + l)-dimensional systems of hydrodynamic 
type which satisfy the Egorov property. Solving these systems by the general- 
ized hodograph method, one can construct exact solutions of the heavenly equation 
parametrized by arbitrary functions of a single variable. We discuss explicit exam- 
ples of hydrodynamic reductions associated with the equations of one-dimensional 
nonlinear elasticity, linearly degenerate systems and the equations of associativity. 
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1 Introduction 



The Plebanski first heavenly equation, 

Vt X yVt Z t Vt X lVt Z y 1, (1) 

governs Kahler potentials of the 4- dimensional metrics 

ds 2 = 2Vt xy dxdy + 2Vt zt dzdt + 2VL xt dxdt + 2Q zy dzdy, 

which are simultaneously Ricci-flat and self-dual |20l l2*T| IT4] . 

This equation was thoroughly investigated during the past three decades using various 
technics. It's integrability via the twistor construction was established in [20] • The 
symmetry properties and conservation laws of the heavenly equation were studied in 
E21 EHl UHl El U] ■ We refer to the review j2] for further references. 

In this paper we propose an alternative approach to the construction of exact solutions 
of the PDE based on its decomposition into a family of commuting (l+l)-dimensional 
systems of hydro dynamic type. This method was successfully applied to the dispersionless 
KP equation in a series of publications [TUl HH H2 UH US] , and subsequently extended 
to the Boyer-Finley equation in [T7J[T31IB]. In principle, it applies to any multidimen- 
sional dispersionless 'integrable' system. Although, written in the form ((TJ, the heavenly 
equation is not quasilinear, one can introduce new variables 

VL xy = a, Vt zt = b, VL xt = p, Vt zy = q, (2) 

transforming ([IJ into the quasilinear form 

a t = p y , a z = q x , b x = p z , b y = q t , ab - pq = 1. (3) 

Hydrodynamic reductions are sought in the form 

a = a(R\. ..,R n ), b = b(R l , R n ), p = p(R l , R n ), q = q^R 1 , R n ), (4) 

where the dependence of the 'Riemann invariants' R 1 , R n on x, y, z, t is governed by a 
triple of commuting hydrodynamic type systems 

R\=v\R)R xl R y = w\R)R x) K z = //(it!) R x ; (5) 

(notice that the number n of Riemann invariants is allowed to be arbitrary!) We recall, see 
|24j . that the requirement of the commutativity of the flows (J5J) is equivalent the following 
restrictions on their characteristic speeds: 



djV 1 djw 1 djfi 1 



V 3 — v % W 1 — W l [I 3 — jl l 



i + h dj = d/d Rj . 



Substituting the ansatz (jlj) into the equations © and using (jSj), one readily arrives at 
a system of PDEs for a,b,p,q,v\w l , // as functions of R 1 ,...,R n , see equations (fT^j) in 
Sect. 2. Although overdetermined, this system is compatible, with the general solution 
depending on 3n arbitrary functions of a single variable. All hydrodynamic reductions of 
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the form (JHJ are automatically semi-Hamiltonian [21], that is, their characteristic speeds 
satisfy the additional relations 



d k f V ■ = 9j ,. , i^j^k, (6) 



(the same constraints holding for w ! and ^*). Relations (JBJ) imply the existence of the 
diagonal metric ^ 9a(dR 1 ) 2 such that 

-d j \n^fg~ u i^j. (7) 



This metric is clearly the same for the whole family of commuting flows (JHJ), 



djV 1 djw 1 djjjL 1 



dj In Jgl 



it plays an important role in the Hamiltonian formulation of the corresponding hierarchy 
[HI El!- The generalized hodograph method proposed by Tsarev [22], provides the general 
solution of the system (JHJ) by the implicit formula 

\\R) = x + v\R) t + w^R) y + /J(R) z, i = l,...,n, (8) 

where X l {R) are the characteristic speeds of the general flow commuting with (jlfij) . that 
is, the general solution of the linear system 

-djhly/jto, I / ./. (9) 



Ultimately, the method of hydrodynamic reductions consists of the following steps: 

- find a, b, p, q,v\w\fi l as functions of R 1 , ...,R n by solving the system (|TT)|): 

— solve the linear system 0; 

- express R 1 , R n as functions of x, y, z, t from the implicit relations (JSJ); 

— reconstruct Q by virtue of (J2J) • 

In Sect. 2 we derive the involutive system of PDEs governing n-component reductions 
of the heavenly equation. We demonstrate, in particular, that all such reductions satisfy 
the Egorov property, that is, the diagonal metric gu defined by (JZJ) possesses a potential: 

ga = dii[). 

Diagonal systems satisfying the Egorov property play a particularly important role in 
applications: examples thereof inlude isentropic gas dynamics, chemical kinetics, equa- 
tions obtained by the Whitham averaging, hydrodynamic reductions of the dispersionless 
KP and Boyer-Finley equations, etc. Egorov's (potential) metrics and the corresponding 
hydrodynamic hierarchies naturally arise in the context of 2D topological field theories 
[7]. Multi-Hamiltonian Egorov's systems were recently classified in |19j . 

One-component reductions (n = 1) of the heavenly equation are described in Sect. 
3. Examples of two-component reductions associated with the equations of nonlinear 
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elasticity are provided in Sect. 4. Some multi-component examples are discussed in Sect. 
5. 

In this paper we concentrate on mathematical aspects of the method. The geometry 
of the examples constructed requires a special investigation. 

Remark 1. More general hydrodynamic reductions can be sought in the form a = 
aiu 1 , ...,u n ), b = b^u 1 , u n ), p = p{u l , ...,u n ), q = qi^u 1 , u n ), where the dependence of 
u 1 , u n on x, y, z, t is governed by a triple of commuting hydrodynamic type systems 

u\ = v){u) u J x , u l y = w){u) u J x , u\ = n){u) u 3 x , 

which are no longer assumed to be strictly hyperbolic or diagonalizable. Reductions of 
this type are presented in Example 3 (elliptic case) and Example 5 (parabolic case). 



2 Properties of n-component reductions 

Substituting a(R), b(R), p(R), q(R) into the equations and using one readily 
arrives at the equations 

diq = fi'dia, dip = —.did, d { b = — -^a; (10) 

w i w % 

the differentiation of ab — pq = 1 implies ad^b + bdia = pdiq + qdip, which, by virtue of 
can be written in the form 

.W^. (11) 

ajjL 1 — q 

Substituting (|TT|) back into (fTU|) . one obtains 

= n l dia, = — — -did, d { b = /i 1 — — -<%a, (12) 
ajjL 1 — q an % — q 



A direct calculation of the compatibility conditions for equations (jl2j) results in the fol- 
lowing over-determined system for a: 

didjd = d jAt g- a + g a 

(13) 

dj/j, l (afi j — q)di<x + diji 1 {aji 1 — q)djd = (// — fi j ) 2 diadja. 
The commutativity conditions, 

djV 1 djW 1 djii 1 ,„ 

- 1 r = — 1 r = 3 . , (14) 

V 3 — V % w J — W l iJ? — fl l 

with v l given by (|llj). imply 

v = v w ^y-.^) 2 _ (is) 

w l (afi J — q) — w^afx 1 — q) 
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Ultimately, combining (fT%Jl and (|15|). we arrive at the following equations for a, q, w l , 



dig = n l did, 




didjd 



diadja(w l +w : >)(fi : > — fj, 1 ) 
w % (afii —q)—wi (afi 1 —q) ' 



(16) 




w l dja(fj, : > — /j 1 ) 2 



w l (afii —q)—wi (afi 1 —q) ' 



i 7^ j; notice that p, b, v l do not explicitly enter these equations. One can show by a direct 
calculation that the system (|T6*j) is in involution, and its general solution depends on 3n 
arbitrary functions of a single variable. For any solution of this system the characteristic 
speeds v l of the first commuting flow are given by the formula (jll)) where the functions 
p(R) and b(R) can be recovered from the equations (|12jl. which are compatible by con- 
struction. Moreover, equations (j!6|) and (|T2*|) are consistent with the algebraic restriction 
ab — pg = 1. Thus, the construction of n-component hydrodynamic reductions of the 
heavenly equation is reduced to the solution of the system (JUJ). 

Remark 2. Introducing a potential ip by the formula diip = dia/w 1 (the compatibility 
is guaranteed by (fTHjl ). one can verify the identity 



implying that reductions (0) satisfy the Egorov property. The corresponding Egorov 
metric is ^2d i ijj(dR' 1 ) 2 . 

Remark 3. We have demonstrated that the system 



which is the third commuting flow in (j^J), satisfies the Egorov property and possesses a 
pair of conservation laws 



such that ab—pg = 1. Conversely, it can be shown that, given a Egorov system with a pair 
of conservation laws satisfying the property ab — pg = 1, one can effectively reconstruct 
the first two commuting flows in (J3J and, thus, obtain a hydrodynamic reduction of (fT|). 
Therefore, the description of hydrodynamic reductions is reduced to the classification of 
Egorov's systems with a pair of conservation laws satisfying ab — pg = 1. Although this 
problem seems to be complicated in general, explicit examples can be constructed by 
studying particular Egorov's systems. As demonstrated in Examples 2 and 3, one and 
the same Egorov's system may possess different pairs of conservation laws satisfying the 
requirement ab — pg = 1. It seems to be an interesting problem to classify such systems 
in general. 

In the two-component case, equations (fT7)|) can be rewritten in the new independent 
variables a and q, assuming the form of a non-homogeneous system for fj}, jj 2 and w 1 , w 2 : 



djfi 



didjip 



0- 



R\ — [i\R) R\ 



a(R) z = q(R) x , p(R) z = b(R) 
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F) 1 _l_ 2a 1 w 1 (fi 1 — fi 2 ) 2 r\ 2 _i_ In 2 ui 2 (^ 1 — im 2 ) 2 

Oafl + H O q fX — „,i( aAt 2_ g )_„,2 (aAt i_ g ), C a /i — ^(a/^-f,)-™ 1 ^ 2 -^)' 

(17) 

a <? «i 1 (a/x 2 — 5)— w 2 {afjb 1 — q) ' a 1 r" 9 «i 2 (a/i 1 — 5)— «) 1 (a^t 2 — g) ' 

Equations ()17)) are a direct analog of the Gibbons-Tsarev system E] governing two- 
component hydrodynamic reductions of the Benney moment equations. 



3 One-component reductions 

One-component reductions of (jHJ) are of the form a = a(R), b = b(R), p = p(R), q = q(R), 
where R(x, y, z, t) solves a triple of hydrodynamic type equations 

R t = v(R) R x , R y = w(R) R x , R z = fi(R) R x , (18) 

which, in the one-component case, are automatically commuting. We recall that the 
general solution of ()18|) is given by the implicit hodograph formula 

f(R) = x + v(R)t + w(R)y + fj,(R)z (19) 

where f(R) is arbitrary. One can readily verify that equations © imply \i = q'/a', v = 
wp' /a', so that the formula (|T9*j) assumes the form 

F(R) = a!{R)x + q\R)z + w(R)(p'(R)t + a'(R)y), (20) 

where F = a'f. Moreover, equations (J3J) also imply that a,b,p,q are subject to the 
constraints 

ab-pq = l, a'b'-p'q' = 0. (21) 

Thus, one-component reductions of the heavenly equation are defined by the formulae (J2J) 
where the functions a(R) , b(R) , p(R) , q(R) satisfy the restrictions (|2*TJl . and R(x,y,t, z) is 
implicitly defined by (|2()j) . These solutions depend on four arbitrary functions of a single 
variable (indeed, F, w and any two of a, b, p, q can be arbitrary). 

Remark. As pointed to us by Y. Nutku, one can slightly 'relax' the scheme described 
above, considering only two hydrodynamic equations 

R t = v{R) Ry, R z = fx(R)R x , (22) 

which, in the one-component situation, are automatically compatible. Indeed, the general 
solution of the system (|2*2"|) is given by the implicit formula 

F{R,t + v(R)y,z + n{R)x) = 0, (23) 

where F is an arbitrary function of three arguments. Then equations Q for a(R), b(R),p(R) 
and q(R) imply 

" = (24) 



along with equations (|21|). Choosing a(R) , b(R) , p(R) , q(R) which satisfy equations (|21|) 
and defining u, p, by the formulae (j23), we reconstruct the solution Q(x,y, z, t) of the 
heavenly equation by virtue of (j2J) where R(x, y, z, t) is defined by the implicit relation 

(J2HJ). 

We emphasize, however, that such 'relaxation' of the general scheme applies to one- 
component case only: one can show that two strictly hyperbolic hydrodynamic type 
systems 

R\ = v\R)R y) Ri = t i i (R)Ri, i>2, 
are compatible if and only if there exists a third system 

R y = w\R) Rl 

which commutes with both. 



4 Examples of two-component reductions 

Example 1. Let us consider the following triple of commuting systems, 

R] = 4e 2R2 -i fil Rl, R 2 = 12e 2R2 -l Rl R 2 X , (25) 

Rl = ie- 4 * 2 ^ 1 Rl, Rl = - 3 - e -* R2 -l Rl Rl (26) 

and 

Rl = -e-^+R 2 ) Rl Rl = e -2(R^R 2 ) R 2. (2?) 
the corresponding metric 

e^+^dR^ + idR 2 ) 2 ] 

clearly satisfies the Egorov property Notice that these systems belong to the hierarchy 
of commuting flows of equations of isentropic gas dynamics with the polytropic constant 
7 = 1. Introducing 

a = — e ^ Rl - 2R2 a = —e-^ 1 - 4 * 2 v =— e ^ Rl+AR2 b = —e~i Rl+2R2 

one can readily verify the identities 0- Therefore, the corresponding solution Q of the 
heavenly equation is given by the formulae 

n xy - 2y _ea , n zt -^3 , n xt _ ^ , n zy - 

notice that these solutions are uniquely characterized by the property Q xy Q zt = 4/3, 
Qxt^zy = 1/3. In fact, one can explicitly calculate all second derivatives of Q in terms of 
R 1 , R 2 as follows: 

O _ 2 2R 1 +2R 2 Q M — ^p^+iR 2 Q _ _J_ p ^R 1 ~2R 2 Q _ 4_fE>l_E>2\ 

3\/3 ' x ~t \/3 ' X U 2\/3 ' xz 

n — 32 p 2 :B 1 +6R 2 o — 4 f 1 R 1 Q R2\ o _ 8 p -2fii +2 R 2 _ 1 2 R^-6R 2 

O _ _^ P -±R 1 -AR 2 q — ?_ P -2R 1 -2R 2 . 

~~ 4^ ' ^ ~~ 375 ' 



the compatibility conditions are identically satisfied by virtue of (|25|) -(|27 |) . Here the 
Riemann invariants y, z, t) and R 2 (x, y, z, t) are implicitly defined by the hodograph 
relations 

X\R) =x + 4e 2R2 -l Rl t + ± e - 4R2 -l Rl y - 2e- 2 ^ Rl+R2 ^ z, 

X 2 (R) =x + 12e 2R2 -l Rl t - le- AR2 ~l Rl y + 2e~^ Rl+R2 ) z, 

where A 1 (i?), \ 2 (R) are the characteristic speeds of commuting flows of systems fl23|) - 
(f2"Tj) satisfying the linear equations 

* A1 i, #4-. m 



A 2 - A 1 ' A 1 - A 2 
The general solution of these equations can be represented in the form 

A 1 = -dx(we-#-#), A 2 = d 2 (we~ Rl - R2 ), 

where the function w solves the linear equation d\&2W = w. 

Example 2. Modifying Example 1, let us consider the following triple of commuting 
systems of hydrodynamic type: 

R\ = e-( Rl+R2 )(cos( J R 2 - R 1 ) - sin(R 2 - R 1 )) i£, 
R 2 = - e -( Rl + R2 )(cos(R 2 - R l ) + sm(R 2 - R 1 )) R 2 X , { ] 



R i = e -(fi 1 +^)( cos ( jR 2 _ R i) + sin(jR2 _ R i^ 

R 2 = e -(« 1 +R 2 )( cos ( jR 2 _ ^1) _ sin(jR 2 _ ^1)) R ^ W 

and 

Rl = -e-^+tf) Rl Rl = e -2(R^R*) R 2 (31) 

Since the z-Qow ()31jl coincides with the z-flow (|27p from Example 1, the corresponding 
Egorov metric will also be the same. Introducing 

a = e Rl+R2 cos(R 2 - R 1 ), q = e^ Rl + R2 ^ sin( J R 2 - R 1 ), 

p = - e Rl+R2 sm(R 2 -R 1 ), b = e^ Rl+R2 ^ cos(i? 2 - R 1 ), 

one can readily verify the identities 0- The corresponding solution Q is given by the 
formulae 

= \e 2 ( Rl+R2 \ Q xt = - e Rl+R2 sin(i? 2 - R 1 ), Q xy = e Rl+R2 cos(R 2 - R 1 ), 
n xz = R 2 -R\ n tt = R 1 + R 2 + sin 2 (R 2 - R 1 ), Q yt = - sin^ 2 - R 1 ) cos(R 2 - R 1 ), 

Q zt = e-iR'+R 2 ) C0S (R^ _ Rl^ Qyy = R 1 +R 2_ sin 2 (jR 2 _ R l^ 

n yz = e -( Rl +« 2 ) sm(R 2 - R 1 ), Vt zz = -\ e -^ Rl+R2 \ 

(32) 



S 



Here the Riemann invariants R 1 (x, y, z, t) and R 2 (x,y,z,t) are implicitly defined by the 
hodograph relations 

A 1 ^) = x + e-( i?1+i?2 )(cos(i? 2 - R 1 ) - sin(R 2 - R 1 )) t+ 
e-( Rl+R2 \cos(R 2 - R 1 ) + sin(R 2 - R 1 )) y - 2e~ 2( - Rl+R ^ z, 

A 2 (#) = x - e-( Rl+R2 ^(cos(R 2 - R 1 ) + sin(R 2 - R 1 )) t+ 
e-( Rl+R2 \cos(R 2 - R 1 ) - sin(R 2 - R 1 )) y + 2e" 2 ( Rl+i?2 ) z, 

where A 1 (i?), \ 2 (R) are the characteristic speeds of commuting flows of systems (|2*9*Jl - 
(j3*T)l which satisfy the linear equations (J28j) from example 1 (indeed, the flows (}25|) . (|26|) . 
IJ27J) . (j2HI), (jHOj) commute with each other). In the new dependent variables 

v = R 2 -R\ f = e Rl+R \ 

equations (f2HJ) - (|3*Tj) take the conservative form 

[f 2 /2] t = \-ip sin v] x , [v) t = [cosv/f] x , (33) 



[f 2 /2] y = [tpcosv] x , [v] y = [smv/(f} x , (34) 

and 

[f 2 /2] z = M„ [v], = [-1/(2^)]*, (35) 

respectively. Written in terms of u = f 2 , equations ()35|) take the form of the nonlinear 
elasticity equation u zz + {l/u) xx = 0. The corresponding equations (J52*|l reduce to 

&xx = If 2 , &xt = - f sin v, Q xy = f cos v, Vt xz = v, 

Q t t = hi f + sin 2 v, Q,yt — sin vcos v, Q zt = cosv/f, fl yy = In ip — sin 2 v, 

Q yz = sinv/f, Q zz = -l/(2ip 2 ). 
Example 3. The elliptic version of equations (|3*3*|) - (j33j) is 

[f 2 /2} t = [ipsmhv} x , [v] t = [coshv/f} x , (36) 



[f 2 /2} y = [tp coshv] x , [v] y = [&inhv/(p] x , (37) 

and 

[ip 2 /2] z = [v] x , [v) z = [1/(2<^ 2 )] X , (38) 

respectively. Written in terms of u = f 2 , equations ()38|) take the form of the nonlinear 
elliptic equation u zz — (l/u) xx = 0. The corresponding equations for Q reduce to 

= f 2 /2, Vt xt = fsinhv, Vt xy = ip cosh-u, Q xz = v, Vttt = sinh 2 t; — Imp, 
Vl yt = sinh v cosh v, VL zt = cosh v/f, VL yy = cosh 2 v + In ip, 
Q yz = sinhv/f, Q zz = l/(2(p 2 ). 
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5 Examples of multi-component reductions 

Example 4. The simplest possibility is to set w l = 1 in equations ()16|). In this case 
equations for a simplify to adidja = 2diddjCi, so that, up to a transformation R 1 — > f t (R i ), 
the function a can be choosen in the form a = l/(-R 1 + ... + R n ). Hence, equations (jlfij) 
become 

a i q = l A t a i a 1 = __!_, (39) 

while v l and p, b are defined by (jlljl and ()12j) . We point out that the choice w l = 1 
implies Q = z,t), £ = x + y, leading to the equation 

%^ 2t - = i, (40) 

which is a reduction of the heavenly equation by a self-dual symmetry [HE]. Therefore, 
equations (jHHJ) (which are manifestly linear) describe hydrodynamic reductions of (f4T)j) . 
Choosing, for instance, 

^ = j e R k - E ° kRk 

where c l are arbitrary constants, one can readily solve the remaining equations (jSHJ)? (jH$ ; 
(fT^j) for g, t>*, p, 6, obtaining 

and 



T, c k R k 



E ' y E ^ fe ' 

P — E ^ ' °-^ K E^ ' 

Thus, we have constructed hydrodynamic reductions of the equation P0|): 

R\ = ( ci E # fe - E c k R k ) R\; 

the associated Egorov metric is 

2 _ (cifi 1 ) 2 + - + (^) 2 

(E^) 2 

The corresponding solution Q(£,z,t) is given by the formulae 

- _J_ O - V R k - (E c k R k )(J2 R k /c k ) 

1 L & — R k ' zt ~ ^ e ^ fc ' 



(41) 



(42) 



(43) 

Here the dependence of R l on £, t, z is defined by the generalized hodograph relations 



q _ T.R h /c k o - E c k R k 
1 _ E « fe ' — E « fe ' 



*w = e + (e ^ - ^ E * + ( c * E ^ - E ° kRk ) 
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where X l {R) are the characteristic speeds of commuting flows of systems (J42j) . that is, 
solutions of the linear system 

W L_ 

\i - A* ^R k ' 

One can show that the general solution of this system is given by the formula 

X i (R)=fl(R i )J2 Rk -J2M Rk ^ 

where fk{R k ) are n arbitrary functions of a single argument. 

Remark. Equations (J42|) can be obtained from the linear systems 

pi ' pi pi J. pi 

.fir — r itr, IX- — C itr, 

t c i z £' 

by a reciprocal transformation 

d£ = (£i? fe ) d£-{J2R k /c k ) dt+ (Y,c k R k ) dz, 

(44) 

t = t, z = z. 

Notice that, by virtue of (|4H|) . the first equation (f4*4"|) implies £ = fi^. This reflects the 
linearizability of the equation (|4T)|) by the contact transformation 

£ = fi^, i = t, z = z, fi = £fi^ — fi, fi| = £, fij = —fit, fii = — fi z , 

transforming (J4"0|) to the linear wave equation 

% + fi fe = 0; 

see, e.g., HE]. 

Example 5. Let us first recall the definition of the equations of associativity (the 
case of 4 primary fields), following [7 . One starts with a function F(t , t 2 , t 3 , t 4 ) whose 
dependence on t l is specified as follows: 

F(t\ t 2 , t 3 , t 4 ) = l -(t l ) 2 t 4 + t W + fit 2 , t\ t 4 ). 

Introducing c a /3 y = dt ^Q^ dtl and rj a p = C\ a p, defining the structure constants = 
rp Cs a j3, and requiring the algebra with these structure constants to be associative, one 
arrives at the following system of third order PDEs for the function f(t 2 ,t 3 ,t 4 ): 

— 2/234 — ^233/223 + ^333/222 = 0, 

_ /244 — /233/224 + /334/222 = 0, 
— /344 — /223/334 + j '224 j '333 = 0, 

(45) 

— 2/234/224 + /244/223 + /344/222 = 0, 
/444 — /I34 + /244/233 — /334/224 + fsAif 223 = 0, 
/333/244 — 2/334/234 + /344/233 = 0; 
11 



subscripts denote partial differentiation with respect to t 2 ,t 3 ,t 4 . Following [7j, one intro- 
duces a triple of hydrodynamic type systems 

t? = <2r& i = 4,tl £ = <&ytl P, 7 = 1,2, 3, 4, (46) 
which commute with each other. Let us take 

/(t 2 ,t 3 ,t 4 ) = ^^-t 4 lnt 3 + t 4 lnt 4 , 

which is a solution of (|45jl. Introducing the notation t 1 = it, t 2 = a, t 3 = p, t 4 = r, one 
can write the flows in the conservative form as follows: 

r t = d x p, Ut = -d x ^r, a t = d x [u+ 2 -fl p t = d x ^ 
d x a, Uy = -d x [± + ^\, a* = + Vy = d x [u+ 2 -fl (47) 



(we point out that these equations are not hyperbolic). One can show by a direct calcu- 
lation that there exists a potential Q whose second derivatives are given by the formulae 

(48) 



^2t — ~2 5 ^j/y — „2 "I" r 5 ^j/2 — _(1 "I" r 2 )> ^22 — ~(1 4" 2 )j 



clearly, Q solves the equation Introducing the variables \l — p/r, v = u + ap/r, the 
general solution of the flows (}4*Tj) can be written in the implicit form 

z = a$ vv - $ M „ + rG v , 
y = $ v + /i[a$ vv - $ M „ + rG v ], 

x — G- n[^r<b w + aG v - + G„], 

where <&(v, fi) and G(v, (i) solve the equations 

&w + = 0, G vv + fi 2 G^ = 0. 

Notice that equations (|4~H|) for the potential fl can be integrated twice, resulting in the 
expression 

n = (£$ w + aG v -^ + G) (v<$> v - $)- 
(a$ vv + rG v - ®n V )(nR + S)- R<S> V + H, 
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where the functions R, S, H are defined by the quadratures 

dR = G^dv - 

dS = (— ~ vG ii)d^ - if^Gfi + vG v )dv, 
dH = [^G, + $„G v ]dv + - *f*]dfi,. 

It is an interesting problem to find other solutions of the associativity equations which 
produce solutions of ((H). 
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